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The induced polarization, P' y , of the neutron in the deuteron photo-disintegration from threshold 
up to 30 MeV is calculated using a variety of different, latest-generation potentials — Argonne vi&, 
Bonn 2000, and Nijmegen I — and a realistic model for the nuclear electromagnetic current operator, 
including one- and two-body terms. The model dependence of the theoretical predictions is found 
to be very small. These predictions are systematically larger in magnitude than the measured P' y 
values, and corroborate the conclusions of an earlier, and much older, study. There is considerable 
scatter in the available experimental data. New and more accurate measurements of the induced 
polarization in the 2 H(7, n) 1 H reaction are needed in order to establish unequivocally whether there 
is a discrepancy between theory and experiment. 

PACS numbers: 21.30.+y,25.20.Lj,24.70.+s,25.10.+s 

I. INTRODUCTION, RESULTS, AND CONCLUSIONS 

The present note deals with the deuteron photo-disintegration and neutron induced polarization in the 2 H(7, ft) 1 !! 
reaction at low energy, from threshold up to about 30 MeV. In this energy range, the photo-disintegration process 
is dominated by the contributions of electric dipole (El) and, to a much less but still significant extent, electric 
quadrupole (£'2) transitions, connecting the deuteron to the np 3 Pj=o,i,2 and 3 Si- 3 Di states, respectively. The 
experimental data are well reproduced by theory, see Refs. 0,0 and Figs'. [DEI The total cross section data are from 

Refs. HSHBSEHminl, ^ angular distribution data at P hoton energy £ 7 =19.8 MeV are from Ref. [3, 
and the data on angular distribution ratios as function of E 1 are from Ref. |l3j |. The calculations shown in these 
figures, various aspects of which are succintly summarized in the following section, are based on a variety of (modern) 
realistic nucleon-nucleon potentials, including the CD-Bonn (BONN) HJ|, Nijmegen I (NIJM-I) 15], and Argonne 
wis (AV18) 0|, as well as on semi- realistic reductions of the AV18 the Argonne vq (AV6) and Argonne v$ 
(AV8) models, constrained to reproduce the binding energy of the deuteron and the isoscalar combination of the 
S- and P-wave phase shifts. In particular, the AV6 ignores spin-orbit interaction components, which are important 
in differentiating among the 3 Po,i,2 channels, and therefore does not provide a good fit to the phase shifts in these 
channels. The model dependence of all theoretical predictions shown in Figs. 2)01 including those corresponding to 
the AV6 and AV8 in Fig. ^ is negligible. 

Calculations of the np radiative capture cross section at thermal neutron energies, based on these same potential 
models, are also found to be in excellent agreement with the measured value, when two-body current contributions are 
taken into account 0,0- The model dependence is again negligible. The np radiative capture up to neutron energies 
of about 100 keV proceeds almost exclusively through the well-known magnetic dipole (Ml) transition connecting 
the 1 So np and deuteron states [l8|. 

On the basis of these facts, one is led to conclude that the Ml and El transition strengths, which the np radia- 
tive capture and deuteron photo-disintegration are selectively sensitive to at low energies, are both consistent with 
experimental data. It is known [l9l | that the neutron induced polarization (Py) in the 2 H(7,n) 1 H reaction originates 
predominantly, in the low-energy regime of interest here, from interference of Ml and El transition terms. Thus, it 
is surprising to find that this observable, measured up to photon energies of 25 MeV, is significantly overestimated, in 
magnitude, by theory, as shown in Figs.0HHl although the data at the (center-of-mass) an gle o f 135° in Fig.[7|seem to 
be consistent with it. The P y angular distribution data at E^—2.75 MeV are from Refs. |20l l2ll |22| . w hile the data 
at (9=90° for £ 7 =6-30 MeV are from Refs. HI El El, and those at 6>=45° and 135° are from Refs~|23i and [HEil, 
respectively. In the figures the results obtained without (IA) and with the inclusion of two-body current contributions 
are displayed separately for the BONN, NIJM-I, and AV18 potential models. The contributions of two-body currents, 
essential if the observed cross section for the np radiative capture is to be correctly predicted, turn out to substantially 
worsen the agreement between the measured and calculated Py in all cases but at 0=135°. The discrepancy between 
theory and experiment is particularly severe for P y at 0=45°. 

That the P' y data are problematic for theory has in fact been known for some time [2^, E^, E3, E3 • Indeed, the main 
motivation for the present study was to re-examine this issue in light of the advances made during the last decade in 
the modeling of both nucleon-nucleon potentials and two-body electromagnetic currents. The corresponding results, 
however, are close to those of Hadjimichael [2?J — reported in Ref. [2j| at 0=90° — and Schmitt et al. [28| and, moreover, 
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show a very small model dependence. There is considerable scatter among the different data sets at 9=90° and in 
the measurements of the P' y angular distribution at -E 7 =2.75 MeV. Clearly, more accurate data on both the energy 
dependence and angular distribution, which could be used to isolate the multipole components, are needed in order 
to resolve this confusing situation, and draw definite conclusions. 



II. CALCULATION 



The relevant matrix element in the photo-disintegration of a deuteron in spin projection m c i initially at rest in the 
laboratory is 



jrn„,m p ;k,m d (P' c l) = P> m »> m P I ^(q) ' j(q) I m d ) (2.1) 

where q is the momentum of the absorbed photon and = ±1, are the spherical components of its polarization 
vector, j(q) is the nuclear electromagnetic current operator, and | q; p, m n , m p )^~^ represents an np scattering state 
with total momentum q and relative momentum p, satisfying incoming wave boundary conditions. The z-axis is taken 
along q, which also defines the spin-quantization axis. In the results of the calculations presented in Sec. [H the np 
state includes all channels up to total angular momentum J=5, the contributions of higher partial waves have been 
found numerically negligible. The methods used to solve for the two-nucleon bound- and scattering-state problems 
as well as the techniques developed for the evaluation of the transition amplitudes above have been described in 
considerable detail in Ref. they will not be discussed further here. 

It is convenient to introduce a second reference frame with axes x', y', and z', in which the relative momentum p 
is along z' with components (sin# cos0, sin# sin0, cos#) with respect to the reference frame defined earlier. The x' and 
y' axes are taken to have directions (cos# cos</>, cos6* sin0, — sm8) and (— sin<^>, cos</>, 0), respectively. A neutron with 
polarization in the +y' direction, as an example, is represented by the state 



l+y'H 1 ' ^ - > ( 2 - 2 ) 

where | ±) denote the spin states with ±1/2 projections along z, i.e. q. The transition amplitude for emission of a 
neutron with polarization in the +y' direction is then obtained from the linear combination 



.•(-) 

J+y' ,m p ;\, 



V2 



'■(-) _i e -i?S-(-) 



(2.3) 



A similar expression holds for emission of a neutron with polarization in the — y' direction. The induced polarization 



P' y is defined as 



= <T+y>(0)-<T-y,(0) 



where the differential cross section is given by 



s -^f - w, t E E i i&^M i 2 • (2.5) 

m d ,X m p 

Here a is the fine structure constant, m is the nucleon mass, and the magnitude p of the relative momentum is fixed 
by energy conservation. The polarization parameters P' x and P' z , proportional to cross-section differences for emission 
of neutrons with polarizations, respectively, in the ±x' and ±z' directions, vanish, as required by parity conservation, 
and this fact has been explicitly verified in the numerical calculations. 

The continuity equation allows one to express the nuclear electromagnetic current entering Eq. (|2.1|l as 



j(q) =j(q)-j(q=o) + i 



H , J dxxp(x) 



(2.6) 
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where H is nuclear Hamiltonian and p(x) is the nuclear charge density operator. This identity ensures that the Siegert 
form is used for the El operator, dominant in the energy regime of interest here. Note that the matrix elements 
J+y' m p -\m d are calculated as discussed in Ref. pj, namely without performing the expansion of j(q) in terms of 
electric and magnetic multipole operators. Of course, the commutator term in Eq. 1|2.1[) reduces to 



i J dxx [H , p(x)] — > iq J dxxp(x) = iqd , (2-7) 

when evaluating the matrix elements. The dominant contribution d NR to the electric dipole operator d, to which the 
(unretarded) El multipole operator is proportional, is simply given by 



d NR = 5^P i (r i -R) , (2.8) 

i 

where R = (ri + r2)/2 is the center-of-mass position vector and Pi is the proton projector. However, there are 
a number of relativistic corrections which have been included in the present study, due to i) the spin-orbit term 
in the single- nucleon charge density operator, ii) the leading two-body contribution to p(x), associated with pion 
exchange, and iii) the center-of-energy correction. Explicit expressions for the associated operators can be found in 
Refs. [23, H3 • In particular, the center-of-energy correction arises because translationally invariant wave functions 
require center-of-energy rather than center-of-mass coordinates. The correct electric dipole operator should be defined 
as 



d = ^P 1 (r,-R CE ) , (2.9) 

i 

where 




(2.10) 



and (for two particles) 



Ei ^m+?L + V -f . (2.11) 
2m 2 



This leads to a correction term to d of the form 



1 CE R — Rce 



8m 2 



l — -(p • P)r + h.c. 



(2.12) 



for transitions between states with zero z-componcnt of the total isospin, P is the pair total momentum. Quanti- 
tatively, however, this as well as the spin-orbit and pion-exchange corrections to d NR have been found to be rather 
small for E 1 up to 30 MeV (see below). It is interesting to note that in the 4 He(d, 7) 6 Li radiative capture, the matrix 
elements of d NR vanish because of isospin selection rules (and the use of translationally invariant wave functions), 
and the relativistic corrections above are resp onsible for the El strength which dominates the cross section for this 
process at energies of 100 keV and below |3fij . 

The two-body currents in the calculations based on the BONN and NIJM-I potentials include the terms associated 
with 7T- and p-meson exchanges, A-excitation, and pnj and ujttj transition mechanisms (see Ref. [l| and references 
therein). The AV18 calculation includes, in addition, the two-body currents associated with its momentum-dependent 
interaction components, as derived in Ref. |2|. It should be emphasized that the AV18 currents are exactly conserved. 

Finally, Fig. [S] is meant to illustrate the sensitivity of the P' y results to various inputs in the calculation. The 
curves labeled AV18, AV8 and AV6 denote the results of calculations based on the AV18, AV8 and AV6 potentials, 
respectively, including one- and two-body currents and the relativistic corrections to the electric dipole operator 
discussed above (the curves labeled BONN, NIJM-I and AV18 in Figs. QHZI are obtained in this same approximation 
scheme, while those labeled BONN IA, NIJM-I IA, and AV18 IA ignore two-body current contributions). As already 
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emphasized, this polarization observable is not very sensitive to the input potential, even in the case of a semi-realistic 
one like the AV6. The curves labeled AV18 IA and AV18 IA w/o RC in El both represent AV18-based results including 
only one-body currents, the difference being that in the AVf8 IA w/o RC in Ef calculation the relativistic corrections 
to dNR are ignored, while in the AV18 IA one they are retained. At the highest energy, the associated contributions 
reduce (in magnitude) the AV18 IA results by about 5% (in the total cross section they decrease the AVI8 IA values 
by less than 1% over the whole E 7 range). The curve labeled AV18 S+P shows the AVI8-based results obtained by 
including one- and two-body currents, but only S- and P-waves in the partial-wave expansion of the final np state. 
As can be seen by comparing the AV18 and AV18 S+P curves, the contributions of the higher partial waves to P' y 
are substantial at the highest photon energies. Lastly, the results denoted as AV18 w/o Ml are obtained with the 
AV18, except that the contribution of the Ml transition due to the state has been neglected. It demonstrates the 
sensitivity of the P' y observable at 0=90° to this component of the amplitude . This sensitivity persists at 6=45° 
and 135°. Note that the total cross section in Figgis reduced by less than 1% for E 1 in the 4-6 MeV range. 
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FIG. 1: (Color online) The deuteron photo-disintegration cross section, calculated with a number of modern nucleon-nucleon 
potentials, is compared to data. Note that the various curves are indistinguishable. 
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FIG. 2: (Color online) The center-of-mass angular distribution of the experimental photo-disintegration cross section, normalized 
to the total cross section, is compared to the results of calculations based on the BONN and AV18 potentials. Note that the 
theoretical curves are indistinguishable. 
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FIG. 3: (Color online) The lab-frame angular distribution ratios, measured in the deuteron photo-disintegration as function of 
photon energy, are compared to the results of calculations based on the BONN and AV18 potentials. Note that the theoretical 
curves are indistinguishable. 
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FIG. 4: (Color online) The center-of-mass angular distribution of the neutron induced polarization measured in the 2 H(7,n) 1 H 
reaction at photon energies of 2.75 MeV is compared to the results of calculations based on a number of latest-generation nucleon- 
nucleon potentials and a realistic model for the nuclear electromagnetic current, including one- and two-body components. Also 
shown are the results obtained by ignoring two-body currents (labeled IA). 
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FIG. 5: (Color online) The neutron induced polarization measured in the 2 H(7,n) 1 H reaction at center-of-mass angle #=90° 
is compared to the results of calculations based on a number of latest-generation nucleon-nucleon potentials and a realistic 
model for the nuclear electromagnetic current, including one- and two-body components. Also shown are the results obtained 
by ignoring two-body currents (labeled IA). 




FIG. 6: (Color online) Same as in Fig. [HI but at center-of-mass angle #=45°. 
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FIG. 7: (Color online) Same as in Fig.Q3 but at center-of-mass angle 9=135° . 
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FIG. 8: (Color online) The neutron induced polarization measured in the 2 H(7,n) 1 H reaction at center-of-mass angle 8=90° is 
compared to results obtained in various approximation schemes (see text for an explanation of the notation). 



